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Abstract 

A new method based on the i?'-operation of the renormalization theory is proposed 
for the numerical calculation of the renormalization constants in the theory of critical 
behaviour. The problem of finding residues of the poles of the Green's functions at e = 0, 
where e = 4 — d, is reduced to the evaluation of multiple UV-finite integrals, which can be 
performed by means of standard integration programs. The method is used to calculate 
the renormalization group functions of the model A of critical dynamics in four-loop 
approximation. Dynamical exponent z of the model A is calculated in the fourth order of 
the e-expansion. 

Introduction 

The method of the renormalization group (RG) nowadays is the main tool of calculation of crit- 
ical exponents in the theory of critical phenomena. The basis of that method is the technique 
of the ultraviolet renormalization of a model. Critical exponents are found from the calcu- 
lated renormalization constants by means of the standard RG rules. Thus calculation of the 
renormalization constants is the primary concern of the method. The most consistent scheme 
of applying the RG method is combining it with the e-expansion (with e = 4 — d, where d is 
spatial dimension) and that is the way it is used throughout the paper. 

Creation of subtle analytical methods of calculations allowed to attain 5-loop accuracy (fifth 
order of the e-expansion) in the theory of static critical phenomena [T] and maximum 3-loop 
accuracy (third order of the e-expansion) in the theory of critical dynamics [2j. These record 
results have held since 1991 and 1984 respectively. Attempts of analytical calculations in higher 
orders met fundamental difficulties [3]. 

For numerical calculation of the renormalization constants it is essential to be able to ex- 
tract residues at e = from the graphs. We developed such technique based on well-known R '- 
operation of the theory of renormalizations. It reduces calculation of renormalization constants 
to the evaluation of multiple UV-convergent integrals which can be performed by standard pro- 
grams for numerical integration. We used the technique for calculation of dynamical exponent 
of the model A in 4-loop approximation (fourth order of the e-expansion). Brief formulation of 
the model and the results obtained are given below. 
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1 Model 



The model A describes dynamical behavior of pure homogeneous Ising system like uniaxial 
magnetic materials. Particularly the critical dynamical exponent z determines the growth of 
the relaxation time near the critical point. This paper presents the results of 4-loop calculation 
of the exponent z with the help of the renormalization group (RG) in the e = 4 — d expansion 
framework. 

The model A in the field-theoretic formulation (see e.g. [3]) can be described by the action 

S(if>, VO = AoVV + $ [-M + Ao(#V - r V - 9oip 3 /Q)] , (1) 

with the n-component order parameter field ?/>(x, t) and the auxiliary field t). At dimension 
d = 4 — e graphs of the perturbation theory are ultraviolet divergent that manifest itself as pole 
at e = 0. These divergences are eliminated by introducing the renormalization constants Z for 
which we use the minimal subtraction scheme (MS). Renormalized action has the form 

S R (i(>, = Z x \$$ + V' [-Z 2 d t ^ + X(Z 3 d 2 ^j - Z^ - Z 5 g^ 3 /6)] . (2) 

It is obtained from (TjQ) by multiplicative renormalization of the parameters Ao = \Z\, tq = rZ T , 
go = gfi e Zg and the fields ip — > Z^ip, ip' — > Z^/ip', with the assumptions Z\ = Z\Z^,, Z 2 = 
Z^'Z^p, Z% = Z^'ZxZ^, Z4 = Z^i Z\Z T Zip, Z<z = Z^i Z\Z g Z^. 

Important extra information about renormalization constants could be obtained from con- 
sideration of the static (simultaneous) Green's functions of the theory (TjQ). It could be shown 
that such functions coincide with renormalized Green's functions of </? 4 theory with static action 

^at ($) = ^,92^2 _ ZiT ^ /2 _ Z5g ^i/24, (3) 

where static renormalization constants Z 3 , Z4 and Z 5 in the MS scheme coincide with the 
respective dynamical ones in §2§ so the renormalization constants of the field ip and parameters 
g and r are equal to the static ones pE]. Moreover the requirement that the static and the 
dynamical theories agree leads to relation 

Z^,iZ\ = Z^ , (4) 

from which it follows 

Zi = Z 2 = Z x x Z 2 , Z^ = Z^p, Z 4 = Z T Z^, Z$ = Z g Z^. (5) 

The only new renormalization constant in addition to the static theory is Z\ that can be 
found from the 1-irreducible functions (ip'ip 1 ) or (ip'd t ip). At present the static renormalization 
constants are known up to 5-loop accuracy (fifth order of the perturbation theory). So in 
the MS scheme one can use coordinate of the fixed point g* (root of /3-function (3{g*) = 0) 
calculated in the static model Q. In the following in place of g we will use for convenience the 
charge u = gS d /(27i) d , where S d = 2vr d / 2 /T(d/2) is the area of the unit sphere in d dimensions. 
The value of the charge at the fixed point u* is given in [3] . With the accuracy required for the 
following calculation we have 

2 34 k 2 2 /-33n 3 + 110n 2 + 1760n + 4544 , . \ 3 4n , n , 

u* = —e + — j{e 2 + — t — 4C(3)A;iA;4 e 3 + 0(e 4 ). (6) 



3fci 81 k\ 21k\ V 5832 
From now on we use the notation: 

k = (3n + 14)/17, h = (n + 8)/9, k 2 = (n 2 + 6n + 20)/27, 

k 3 = (n + 2)/3, k 4 = (5n + 22)/27, (7) 

all ki are normalized so that ki — 1 at n = 1. 
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2 Calculation of the dynamical critical exponent 

The dynamical critical exponent z is expressed by the relation z = 2 + 7* — 77 by means of the 
Fisher critical exponent 77 and the value 7* = j(u*) of RG-function j(u) = P(u)d u \n Z 2 at the 
fixed point u* [I]. Traditionally z is written in the form z = 2 + Rr) with R = 7*/^ — 1. The 
quantity R is convenient because the first two terms of its e-expansion being independent of n 
(the number of field components). Indeed the leading terms of the ^-expansion of the Fisher 
exponent are determined by the series jji] 

u* 2 

V = h— [1 + axhu* + (a 2 k 2 + a 3 k 3 + a 4 k 4 )u* 2 ] + 0(w* 5 ), (8) 

where 

ai = -3/8, a 2 = -15/64, a 3 = -5/32, a 4 = 45/32, (9) 
while the analogous expansion for the 7* is given by 

u* 2 

7* = hh— [I + c x k lU * + {c 2 k 2 + c 3 k 3 + c 4 k 4 )u* 2 } + 0{u* 5 ), (10) 



with some constants q, h. The quantities ([8]) and (flOl) depend on n only via ki. Calculating 
the quotient 7*/?? from these relations and using ([6]) one finds 



7* f 2 

— = h{ i + -( Cl -ai)e + 



4 4 nh + k 

— ai ai -ci H — Co — ao + tt: 

9 v ' 3 V ; (n + 8) 2 



e 2 + 0(e 3 ) , (11 



where 



176 88 
b = 28(ci - ai) — (c 2 - a 2 ) + 24(c 3 - a 3 ) + —(c 4 - a 4 ), 

40 20 
61 = 6(ci - ai) - y (c 2 - a 2 ) + 12(c 3 - a 3 ) + y(c 4 - a 4 ). (12) 

It is seen from these relations that first two terms of the e-expansion of R indeed do not depend 
on n. Evaluation of the renormalization constant Z 2 in 4-loop approximation gives the following 
values of the coefficients in ffTUl): 



h = 1.72609, ci = -0.4939(4), c 2 = -0.2512(0), c 3 = -0.1699(4), c 4 = 1.806(3). (13) 

The first calculation of h was performed in [5] (one 2- loop graph): 

/i = 61n(4/3). (14) 

The constant C\ was calculated in |6J and corrected in [2] (three 3- loop graphs): 

tt 2 /8-F(1/4) 3 13 21 

ci = ' , \' -- + — n4 n3, 15 

ln(4/3) 4 8 8 V ; 

where F{x) = dt\nt/(t — 1) is dilogarithm (Spence's function). The values of the constants 
h and c\ obtained in ( FTBl agree with the analytic expressions from f Tl4"l) and ( |T5l) . 
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In regard to the constants we pioneer in calculating of, namely, c 2 , C3 and C4 (twenty five 
4- loop graphs). It turns out that one can verify the value of c 2 using the asymptotic expression 
for the constant R in the limit n — > 00 [5]: 

fl= 4f (4-^(l- e /2) 1 

Calculating of the ^-expansion in the rhs of ( fl6l) gives 

i?= [61n(4/3) - 1](1 - 0.188483£- 0.0999529£ 2 ) + 0(e 3 ), n -> 00. (17) 

On the other hand from (Hip we obtain 

7* f 2 4 1 

i? = — -1 = hll + -(ci -ai)e + -K (ai - c x ) + 3(c 2 - a 2 )]e 2 j -l + 0(e 3 ), n -> 00. (18) 

After the substitution of the values from and (I13p into (1181) the agreement between the 
coefficients of the e-expansions in (fTTI) and (TTH|) is established. 

The final expression of the constant i? is obtained by the substitution of the quantities (TTcfl) 
to (T2D and (TUD: 

4.78(6) ^21.5(4) 

(n + 8) 2 V ; 



R=[6 ln(4/3) - 1] (l - 0.1884(9) e + 



e^+Oie 6 ). (19) 



The final expression for the dynamical critical exponent z = 2 + Rr] for n — 1 results from (j4j), 
P) and (P2j): 

z| = 2 + 0.0134462 e 2 + 0.011036(2) e 3 - 0.00558(5) e 4 + 0(e 5 ). (20) 



3 Discussion 

The method that we used for extracting the residues of graphs at e = and for the numerical 
calculating of the renormalization constants, turned out to be highly effective for the RG- 
analysis of the dynamical model A in the 4-loop approximation. Correctness of the obtained 
results was confirmed by the cancellation of the poles at e = in the RG function 7(2/) = 
P(u)d u In Z 2 and by comparison with the results known in two and three loops and with the 
leading order of the 1/n-expansion. The main difficulty in the approach we used was achieving of 
high precision in numerical calculations of multiple integrals (in our case maximum multiplicity 
equals to nine and the achieved precision is of order 0.1%.) The proposed method can be used 
for multiloop calculations of the renormalization constants for both dynamical and statical 
problems. In particular, we are planning 6- loop calculation in the static model (j3J). 

It is obvious from (1201) that the 4-loop correction to the critical exponent z is relatively 
small for the real value of e — 1 (i.e. d — 3) we obtain z = 2.0188(9). 
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